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Abstract 

We give a systematic treatment of lattice Green functions (LGF) on the d-dimensional 
diamond, simple cubic, body-centred cubic and face-centred cubic lattices for arbi- 
trary dimensionality d > 2 for the first three lattices, and for 2 < d < 5 for the 
hyper-fcc lattice. We show that there is a close connection between the LGF of the 
d-dimensional hypercubic lattice and that of the (d — l)-dimensional diamond lattice. 
We give constant-term formulations of LGFs for all lattices and dimensions. Through 
a still under-developed connection with Mahler measures, we point out an unexpected 
connection between the coefficients of the s.c, b.c.c. and diamond LGFs and some 
Ramanujan-type formulae for \/n. 

1 Introduction 

Lattice Green functions arise in numerous problems in condensed matter physics, such 
as lattice vibration problems, luminescence, diffusion in solids and the dynamics of 
spin waves [5j . They are also central to the theory of random walks on a lattice [3D] , 
and to the calculation of the effective resistance of resistor networks [T7] . 

In the framework of random walks, for a translationally invariant walk on a d- 
dimensional periodic Bravais lattice, a natural question to ask is the probability that 
a walker starting at the origin of a lattice will be at position I after n steps. The 
probability generating function is known as the Lattice Green Function (LGF). It is 

The coefficient of z n , denoted [z n ]P(l;z), is the probability that a walker starting 
at the origin will be at position I after n steps. Here A(fc) is the structure function 
of the lattice, and is given by the discrete Fourier transform of the individual step 
probabilities. Note that the structure function is not unique, as a given lattice may 
usually be represented by more than one set of basis vectors. We give a simple example 
of this in eqn © below. 

For example, for the d-dimensional hypercubic lattice, the structure function is 
X(k) — ^(cos ki + cos hi H cos ~k&). The probability of return to the origin is 

1-1/P(0;1). 
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In two dimensions we have 

1 r r dk 1 dk 2 



(2tt) 2 J_ n J_ n l-\(k) 

Since P(0; 1) diverges for two-dimensional lattices, this leads to the well-known result 
that the probability of return to the origin by a random walker in two dimensions is 
certain. 

Of broader interest are LGFs defined by: 

,~ , i r r d d % 

The coefficient [z"]P(0; z) is the probability that a walker starting at the origin returns 
there after n steps. For a regular lattice of co-ordination number q, an equivalent 
statement is that q n [z n ]P(0; z) is the number of returns to the origin of an n-step 
random walker (and is of course an integer). For the regular two-dimensional lattices 
the structure functions ar45 

>>(k) honeycomb = + 4 COS 2 k x + 4cOsfci COSfc 2 ). 
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K k )square = -(cosfci +COSfc 2 )- 



K^triang = ^(cosfci + COS fc 2 + COS (k X + fc 2 )). 



The corresponding LGFs are: 



P(0; z) honev = 6v ^ =K(fc) (3) 

V )honey y ) y , 



where 

Az 2 

K = 



(3-z) v /z(3-z)(I + z) 

and K(fc) is the complete elliptic integral of the first kind. For readers who prefer 
hypergeometric representations, recall that 

K(fc) = | 2 F 1 (i \;l;k). 
For the square lattice, the result is remarkably simple, 

P(0;z) sq = -K(z), (4) 
while for the triangular lattice it is a bit more complicated: 

P(0; z) tri = -^K(fc') (5) 



1 Because the honeycomb lattice has two types of site, the expansion parameter z in eqn ((2)1 should be 



replaced by z 2 . 
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where c = (a + l)(b — 1), and 



a = - + 1- \/3+-, and b=- + l + \3+- 



and 

l2(b-a) 



k' 

c 

For the square lattice, we can also use the equivalent structure function 

Hk) square = cos/ci cosfc 2 . 

as the square lattice can be considered either as the two dimensional hyper-cubic lattice, 
which gives the first form of the structure function (above), or as the two-dimensional 
hyper-body-centred cubic lattice, giving the second form. While the integrands are 
clearly different, the integrals are equal, demonstrating the point made above that 
structure functions for a given lattice are not unique. We have 

dk 1 dk 2 r dk x dk 2 

z(coski cos k 2 ) ' 

(6) 

Similarly, for the honeycomb lattice we can exploit the duality with the triangular 
lattice and write the structure function as 



, , r r dk Y dk 2 r r 

4tt 2 P 0; z ) sqU are - / / z ~i Hi TT = / / T 

J-nJ-n 1 - § (COsh +COsfc 2 ) J _ n J _ n 1 



2/1 \ 

^(k)honey "3(2^ ^(^tri J ■ 

dk\dk 2 
1 + 

1 r r dk x dk 2 



It follows that 

P(0;z) honey = (0^)2 J_ v 1 _ ^ ( i + | [c os fcf+cTsfe + cosfa + fc 2 )]) (7) 



(2-7r) 2 J_ n J_ v 1-^.(1 + 4 cos 2 fci + 4cosfci cosfc 2 ) ' 

It is also instructive to consider the coefficients in the expansion of the LGFs, as 
these immediately give the number of returns to the origin after a given number of 
steps. We have: 

P(0;z) = J2*n£) n (8) 

n>0 q 

where q is the co-ordination number of the lattice. Thus q = 3 for the honeycomb 
lattice, q = 4 for the square lattice and q = 6 for the triangular lattice. For the 
honeycomb lattice, 

2 n / \ 2 /ri ■ 

n\ I 2 j 



. j\k\l\ I *r~L \j J V i 

j+k+l=n VJ 7 J=0 vy/ w 



— £ (jSh) =£(") (7)- < 9 > 

For the square lattice, 



a 2 „ = 



3 



For the triangular lattice, 

a n = ( . J (— 3) n ~^bj, where bj — a% 3 (honey comb). (11) 
3=0 ^ 

The result for the triangular lattice appears to be new and has been derived by 
exploiting the connection between the triangular and honeycomb structure functions, 
following a similar calculation in |13j for the diamond-f.c.c. pair. 



1.1 Watson integrals; z — 1. 

Setting z = 1, P(0; 1) gives the famous Watson integrals for the d — 3 case. These were 
first encountered by van Peyp^l, a student of Kramers, who solved the b.c.c. case, but 
Watson solved the more difficult s.c and f.c.c cases, and this family of integrals now 
bears his name. The structure functions are: 

A(fc)<jiom = j(l + (cos fci cos k 2 + cosfc 2 cos/c 3 + cos ki cosfc 3 )). 

X(k)sc — o( C0S ^1 + C0S ^2 + cos ^3). 
o 

A(fc)bcc = (cos ki COS k 2 COS A3). 

A(fc) / cc = -(cos fex cos k 2 + cos &2 cos fc 3 + COS fex cos k 3 ). 



The Watson integrals are given by: 

1 r r r dk x dk 2 dk 3 



P(0;1) 



{2tt) 3 J_ W J_ W J_ W l-A(jb) 
and have been evaluated to yield: 

P(0; l) dwmo „d = -^[r(l/3)] 6 = ^P(0; l) /cc w 1.79288 
^(0; 1). = i(V3 - l)[r(l/24)r(ll/24)] 2 a 1.516386 
P & = ^ii73^4 [ r (V4)] 4 « 1-3932039 

^(0;l)/ cc = ^3 [r(l/3)] 6 « 1.344661 

The result for the simple-cubic case is a saga in itself. Watson gave the result in 
terms of the square of a complete elliptic integral with an algebraic argument in 1939. 
Some 40 years later Glasser and Zucker [26] showed that Watson's result could be 
simplified to 

P(0;i)sc = ^[r(i/24)r( 5 /24)r(7/24)r(ii/24)] 2 , 

which J Borwein and Zucker simplified to the quoted result some 15 years later in |11) . 



2 An even earlier mention is by R Courant in 1928 in the simple-cubic case. 
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1.2 The general case for arbitrary z. 

More generally, the LGFs for the three-dimensional lattices are also known. For the 
simple cubic lattice one has: 



P(0; z) 



1 



7T pTT pTT 



dk\ dk 2 dk% 



l 71 ") 3 Jo Jo Jo 1 - f (cosfci + cosfc 2 + cosfc 3 ) 
Joyce [33] showed that this could be expressed as 



P(0;z) 



1 - 9^ 4 



where 



(l-£) 3 (l + 3£) 
16£ 3 



: K(A;i) 



(l-£) 3 (l + 3£)' 



with 



e = (i + Vi-z 2 y 1/2 (i - vi-^ 2 /9) 1/2 - 

For the body-centred cubic lattice one has: 



P(6;z) 



7T pTT p 7T 



dki dk 2 dk 3 



i 77 ) 3 Jo Jo Jo 1 — 2: (cos k\ COS k 2 COS ks) 
Maradudin et al. [37] showed that this could be expressed as 



P(0; z) = 



-K(fca) 

7T 



where 



1 1 



2 2 

For the face-centred cubic lattice one has: 

p°7T pTT /*7T 

P(0-z 



JO JO 



cZfci dfc 2 <^3 



1 



|(CiC 2 "I- C1C3 + C2C3) 



where c, = cos fcj . Joyce [33] ) showed that this could be expressed as 

n 2 



nu ' z) ( i-e) 3 (i + 30 



;K(fc 3 ) 



where 



16£ 3 



(i-0 3 (i + 30' 



and 



e = (1 + + s/T+Tz). 

Finally, for the diamond lattice one has: 



P(0;z) = 



7T pTT pTT 



dk\ dk 2 dk% 



O) 3 Jo Jo Jo 1 - %(1 + cic 2 + cic 3 + C2C3) 
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where 



which Joyce [3T] pointed out gives 

F(0;z) = 4k(^)K(L); (12) 

fc| = I ± 1*2(4 - z 2 )(V2) _ 1(2 _ z 2 )(1 _ JS 2j(l/2)_ 

The derivation of these results and some history and erudite discussion can be 
conveniently found in the book by Hughes [30]. Note too that as the LGFs can be 
expressed as the square of a complete elliptic integral, they can also be expressed as 
the square of a 2-F1 hypergeometric function. Then by a Clausen type formula they 
can be expressed as 3-F2 hypergeometric functions, as shown by Joyce |32j . We give a 
slightly simpler representation for the diamond and f.c.c. lattices in Section 3.1 

As discussed above for the two dimensional lattices, it is informative to consider 
the coefficients in the expansion of the LGFs, as these give the number of returns to 
the origin with a given number of steps. We have: 

P(0;z) = E a »0" ( 13 ) 

n>0 9 

where q is the co-ordination number of the lattice. Thus q — 4 for the diamond lattice, 
q = 6 for the simple-cubic lattice, q = 8 for the body-centred cubic lattice and q = 12 
for the face-centred cubic lattice. For the diamond lattice, 



l 2n 



j+k+l+m—n v j— 

For the simple cubic lattice, 



E (jW.Uml) (/)( „_/)" (14) 



= /2n\ y, fj±_\ _ (2n\ ^ (n\ (2j 

^ 71 ' j+k+l= 

For the body-centred cubic lattice, 



ogT' = ( 2 ;) (16) 
while for the face-centred cubic lattice, it was recently pointed out in |13j that 

^ cc) =E(-)(- 4 r j ^ (17) 

where bj — a^j , given in eqn (fT4|) . 

The results presented to date are generally well-known, apart from (|11|) which is 
new as far as we know, and the analogous result (|17p which has only recently been 
obtained. In the next sub-section we discuss the form common to 3-dimensional LGFs, 
and in the following sections we discuss higher dimensional LGFs which then allows us 
to consider them in a more general setting, with dimensionality as a parameter. 
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1.3 Properties of LGFs for d < 4. 

In summary, we see that the two-dimensional LGFs for the standard lattices can be ex- 
pressed in terms of the complete elliptic integral of the first kind. For three-dimensional 
lattices, the LGF can be expressed as the product of an algebraic expression and two 
elliptic integrals of the first kind, with moduli that are algebraic functions of the ex- 
pansion parameter z. In [35] Glasser has shown that this is a consequence of the fact 
that the LGFs for three dimensional lattices can be reduced to the following form: 

f K(M) du = -?-[V(l-o)(l + ft) - V(l + a)(l-&)]K(fc+)K(fc_) 
J a y/(b-u)(a-u) °-a 

where k± is an algebraic expression of a and b. 

Another way of looking at the LGFs of the three-dimensional lattices is in terms of 
the underlying ODEs. For example, for the s.c. lattice one has [3T] 

4x?(x-l)(x-9)P"' {x)+12x(2x 2 -15x+9)p" (x)+3(9x 2 -Ux+12)P (x)+3(x-2)P(x) = 

(18) 

where x — z 2 . This ODE can be written in the form eqn (|52[) required by Appel's 
reduction [7] which gives the solution of certain third-order ODEs in terms of the 
squares of the solution of an associated second order ODE, notably ([53]). It turns out 
that all the 3d LGFs satisfy ODEs with this very special property. The corresponding 
ODEs for the other lattices are given in Joyce [31]. In the case of the s.c. lattice, the 
associated second order ODE is 



V 0) + 



1 1 1 



x 2(x-l) 2(x-9) 



V 0) 



3(>-4) 



16jc(x- l)(x-9) 



2/0) =0. (19) 



In this form, the ODE is readily recognised as the Heun equation, and indeed the 
solution of psp analytic around the origin is 

P sc (z) = [P(9, -3/4; 1/4, 3/4, 1, 1/2; z 2 )] 2 (20) 

where F is a Heun function. 

The LGFs of the other three-dimensional lattices (b.c.c, f.c.c, diamond) are also 
given by third order ODEs which satisfy the Appell condition. The ODEs can be 
found in Joyce [34 . For the b.c.c. lattice the reduction to a second order ODE results 
in the simpler solution in terms of the hypergeometric function, 

Pbcc(z) = [ 2 Px (1/4, -1/4; l;z 2 )] 2 . (21) 

For the f.c.c. lattice one again has a solution in terms of the square of a Heun 
function, 

P /cc (z) = [P(-3,0;l/2,l,l,l;z)] 2 , (22) 
while for the diamond lattice the situation is similar, 

P d {z) = [P(4, -1/2; 1/2, 1/2, 1, 1/2; z 2 )] 2 . (23) 
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As an aside we remark that the Heun functions above can all be expressed as a 
product of an algebraic function and a hypergeometric function, with a complicated 
algebraic argument. So for example the diamond lattice solution (|23|) can be written 



P d (z) = [y/l - z 2 - yj\ - z 2 ][ 2j F\(1/2, 1/2; 1; k d )f 

where 



with similar results pertaining for the s.c. and fee. lattices. This hypergeometric 
function is immediately recognizable as — K(fci). 

The value of the LGF at z = 1 gives the Watson integrals, and they can all be 
expressed as powers of Gamma functions with rational arguments, times a numerical 
factor involving powers of tt. In two dimensions the analogous integrals are divergent. 

The coefficients in the series expansion of the LGF can be expressed as either 
powers of a binomial coefficient in the simplest case, or as sums of products of binomial 
coefficients, or as double sums of products of binomial coefficients. 

Most attention has been devoted to the calculation of the Watson integrals P(0; 1), 
but in some circumstances, such as the evaluation of effective resistance in networks of 
resistors, it is useful to consider the LGF P(l; 1). This has been done, by Glasser and 
Boersma [23], who showed that 

P(h 1) = A ■ P(0; 1) + - — + C, 

V ' V ' ^ 2 P(0;1) 

where A, B, C are lattice-dependent rational numbers, which depend on I, and satisfy 
simple recursions. For a discussion of resistor networks using LGFs, see also Cserti in 

D3I- 

Other generalisations include the anisotropised version of the LGF, in which the 
structure function for, say, the s.c. lattice is generalised to 

A(fc) = (acosfci + cosfc2 + coskz). 

2 + a 

The corresponding LGF is evaluated both at and at other points in a tour de force by 
Delves and Joyce |18j . Other anisotropic structure functions, or linear combinations 
of structure functions from different lattices, have also been variously treated, see for 
example Joyce, Delves and Zucker in [35]. In all cases the results can be expressed as 
a product of two elliptic integrals. Some insight into why this is so can be found in 
Glasser's calculation in [22] , discussed above. 



1.4 Connections with number theory. 

In a recent paper, Rogers |41| obtained new Ramanujan-type formulae for 1/tt by 
considering the logarithmic Mahler measure of certain 3-variable Laurent polynomials^. 
Mahler measures arc defined and discussed in Section 3.1 below. Theorem 3.1 of Rogers 

3 In fact the second identity in Rogers's theorem, giving eqn H25|) . appeared previously in [16] . 
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[41] gives a number of new formulae for 1/tt in terms of sums involving products 
of binomial coefficients. From equations (fl4l [15]) above, these can be recognised as 
coefficients of the LGF of the diamond and simple cubic lattices. Accordingly, we have, 
as corollaries, the following surprising formulae for 1/tt in terms of LGF coefficients. 

71 — 

8%/3 _ ^ (5n + l) ( d ) , . 



n=0 



5^ 



n=0 




5^(6n + 3-V3) -^-^ 4 d l (26) 



2 ( 64 + 29 ^ = £>20» + 159 - 48V3) ( 80 ^ 139< 

n=0 \ / 



(27) 



One of Ramanujan's 17 formulae for 1/tt [30] can also be expressed in terms of the 
number of returns on the b.c.c. lattice. 

4 _ ^ (6n + l) (6ec) 

--2.^^^ • (28) 

71=0 

There is also the remarkable Ramanujan type formula due to the Borwein brothers [TUJ 
which allows one to calculate any binary digit of 1/tt without calculating the earlier 
part of the binary expansion. This can be written as: 

16 = g(42, + 5) 4&) 

tt ^ 4096" 2n y ' 

n—0 

Note that the general form of a Ramanujan type series for 1/tt is as follows [45] : 
af(z Q )+pef(z ) = -, B = z^- (30) 

TT dz 

where a, /3 and Zq are algebraic numbers and f(z) is the analytic solution around 
the origin of a particularly "nice" third order Fuchsian linear ODE. As pointed out 
by Zudilin [45], this ODE must be the (symmetric) square of a second order ODE, 
a situation first discussed by Appell [7] more than a century ago, and explained in 
equations ([52l [53]) in Appendix A below. As noted in the preceding section, the three 
dimensional LGFs for all four common 3-d lattices satisfy the Appell property, so in 
that sense it is not surprising that they are candidates for a Ramanujan type formula. 
What is more surprising is that they actually do occur in this way 
For example, the formula given by eqn (|2"7) satisifes (|3"01) with 



a = 
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To end this section we point out that there exists another family of formulae for 1/tt 2 , 
called the Ramanaujan-Guillera formulae. Zudilin |45j has shown how these may be 
constructed from Ramanujan type formulae for l/ir by taking the fourth power of the 
solution of the associated second-order ODE (1551) that appears in Appell's reduction. 
In this way we could construct known formulae for l/ir 2 from the alternative route via 
the LGF ODEs. 

2 Lattice Green Functions for d > 3 and Calabi-Yau 
differential equations 

Calculating the LGF for d > 4 is a simple exercise in the case of the hyper-bec lattice, 
and rather more difficult in the case of the hyper-f.c.c lattice. In the case of the four- 
dimensional lattices, a unifying feature is that the LGFs all satisfy a 4th order ODE 
of Calabi-Yau (C-Y) type. The definition of C-Y ODEs and their properties is given 
in Appendix B. 

Calabi-Yau equations are ODEs that have the formal properties of operators that 
appear as Picard-Fuchs operators. More formally, consider a family of n-dimensional 
CY manifoldijf) X(z) parameterized by z G C, where C is an open set on a compact 
Riemann surface. By definition, X(z) has a unique holomorphic n-form w(z) that 
vanishes nowhere. The periods f2j(z) = J r oj(z) are the integrals of oj(z) over the 
n-cycles I\ that form a basis for the n-th homology of X(z). As we vary z on C, each 
period fli(z) will vary in such a way that it satisfies an ODE in z. This is the Picard- 
Fuchs equation for the period of X(z). For further details, see for example |39j . 
CY ODEs comprise a class of ODEs that are pivotal in string theory. 

Here we consider only 4 th order ODEs, (corresponding to the case of Calabi-Yau 
threefoldsjl In [21] the authors define C-Y ODE's as (4th order) ODEs satisfying 
five conditions. The list of known equations is given in [3], and in an updated form at 
http : / / enriques . mathematik. uni-mainz . de / C Yequations / 

Consider ODE's of the form 

y (s) + a s _i(z)y( s ~ 1 ) + • • • + ai {z)y' + a (z)y(z) = 0, 
where {a^} are meromorphic fns. of z. If z — is a regular singular point, we can write 

a s -j(z) = z~ j a s -j(z) j = l,..., 8, 
where a s _j(z) are analytic at z = 0. Then the roots of the indicial equation 
A(A - 1) •• • (A - s + 1) + a s _i(0)A(A - 1) • • • (A - s + 2)+ 

+ ---+5i(0)A + a (0) = 

4 A Calabi-Yau manifold of dimension n is usually defined as a compact n-dimensional Kahler manifold 
satisfying certain technical conditions. A one-dimensional Calabi-Yau manifold is a complex elliptic curve, 
and is algebraic. 

All ODEs considered in this article are linear. To save repetition we delete that adjective, but it should 
be understood to be implied 
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determine the exponents of the ODE at the origin. 

The first condition is that the ODE must have Maximal Unipotent Monodromy 
(MUM). An ODE has MUM if all the exponents at are zero. 

Consider a 4 th order ODE which is MUM: 

+ a 3 (z)y (3) + a 2 (z)y" + ai (z)y' + a (z)y = 0. 

It has four solutions, yo, yi, 2/2, J/3- 
Being MUM implies that: 



2/1 



1 + J2 a nz n ; (31) 

n>l 

2/o log z + ^ b n z n ; 



n>l 



!)2 = -;/„ loj;" : [ y I/,,-:" | log z + ^ c„z 

n>l 



2/3 

The second condition is 



iy log 2 z + ^ 6„z" j 

log 3 z + i ( b ^ n lo s 2 z + ( E c « z " lo § z + E d » z "; 



ai = -a 2 a 3 - -a 3 + a 2 - -a 3 a 3 - -a 3 . (32) 

The third condition is that the roots of the indicial equation at z — oo, Ai < A 2 < 
A 3 < A4 are positive rational numbers satisfying Ai + A4 = A 2 + A3. 

The fourth condition is that the power series solution y has integral coefficients. 

The fifth condition is that the genus zero instanton numbers (defined below) are 
integral, up to a multiple of a fixed positive integer. 

Define q = exp(j/i/j/o) = J2 n >i tn,z n ; the inverse function z = z{q) = ^ u n q n is the 
mirror map in C-Y language. Then the Yukawa coupling K(q) is given by 

d YfvA 1 , ^ fc3jv fcg fe 

%) u + §w 

iVfc are called instanton numbers, and two C-Y equations are considered equivalent if 
they have the same instanton numbers. 

The combinatorial or probabilistic significance of the coefficients of K(q) or of the 
instanton numbers is not known, except insofar as they are of "geometric origin" , as 
discussed in Appendix B. 

Almkvist et al. [3] have catalogued a large number of 4 th order ODEs. They, and 
we, use the operator 

a d 

az 



and write the 4 th order ODE as 

=0, 



11 



where 



V = 6 4 + zPxiB) + z 2 P 2 {6) ... + z k P k {9) 

where P;, / = 1 . . . k are polynomials of degree 4 in 9. 
Thus ODEs of 1st degree take the form 

[9 4 + zP 1 (9)}f(z) = 0. 



They are ODEs of generalised hypergeometric functions. Almkvist et al [3] found 
exactly 14 such ODEs. All can be solved. Their solutions have coefficients expressed 
as finite sums of products of binomial coefficients. We now consider the C-Y ODEs 
that arise in the 4-dimensional generalisations of the four common lattices discussed 
above. 

2.1 Hyper body-centred cubic lattice 

For the 4-dimensional hyper-bec lattice, the LGF is 



Expanding 1/(1 — Xz) as a power series in z, which is absolutely convergent for 
\z\ < A, we have 



This admits to no further simplification, or special values at z = 1, though the 
hypergeometric function is rapidly evaluable, giving 

P(0; 1) fa 1.1186363871641870683496192575256409167948575515294. 

The LGF satisfies a first-degree, 4th order C-Y ODE which is number 3 on the list 
of Almkvist et al. 

With 6 = z£, the LGF satisfies VP(0; 16 z) = 0, where 



Because of the simple structure, we can also easily treat higher dimensions. For 
arbitrary d we have 
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_ ^ / 2 n\ d /_z\2» 
~^{n) \2 d J 

n=0 v 7 

which satisfies a d th order Fuchsian ODE, 

(6 n - z2 n (l + 26>)") P d (6; 2 d z) = 0. (33) 

For d = 5 this equation is of the form required by the Yifan Yang pullback [2], 
which means that the solutions of the 5 th order equation can be expressed in terms 
of the solutions of a 4 th order equation. This is made more precise in Appendix A. 
For this particular ODE, given by eqn (|3"3"|) above, at order 5 (corresponding to the 
5-dimensional hyper bcc lattice), the details of the YY-pullback are precisely worked 
through in Section 7 of Zudilin's article [35]. The pullback results in equation number 
204 in 4 J. A similar result is true for the 5-dimensional hyper cubic (|37p and hyper 
diamond (|39|) lattice LGFs, as discussed below. That is to say, the 5th order ODEs 
permit a YY-pullback to 4th order ODEs which are in the cited list. 



2.2 Hyper cubic lattice 

For the 4-dimensional hyper-sc lattice, the LGF is 

P{ 6., z)= if\..r ******** (34) 

n 4 Jo Jo l-f(ci+c 2 + C3 + c 4 )' 

where Cj = cos fe, . 

Following the treatment in [23] , we use the identity j = J °° exp(—Xt)dt, to write 
the LGF as 



= f «-''o 4 (f)=x>»* 2 ". 

"'O ^ 7 n=0 

since I (z) = ± J* e zcose d6. 

Clearly the use of this identity reduces the d-dimensional hypercubic LGF to a 
one-dimensional integral, involving the d th power of the Bessel function Io(zt/d). 

In this way, Glasser and Guttmann obtained a series expansion which was used to 
identify the underlying 4th order ODE of degree 2. It turns out to be equation number 
16 on the list of Almkvist et al. With = z£, the LGF satisfies VP(0; 8z) = 0, where 

V = A - 4z(26 + l) 2 (56 2 + 56 + 2) + 
+ 2 8 z 2 (6 + l) 2 (26 + l)(26 + 3). 
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The solution is given in [1] as 




The first equality follows from a result [25] of Glasser and Montaldi, who related the 
hypercubic LGF to the generating function of the squares of multinomial coefficients, 
which had been previously given by Guttmann and Prellberg [29] . 

The second equality is in Almkvist [3J, who gives a further 15 distinct expressions 
for a n , all involving single or double sums of products of binomial coefficients. 

Almkvist (private communication) also pointed out that the Yukawa coupling in 
this case is 

K(q) = 1 + 4q + 164,? 2 + 5800g 3 + 196772 g 4 + • • • , 
and hence the instanton numbers are: 

{3N k } = {12, 60, 644, 9216, 157536, 3083604 • • • } 

From the work of Glasser and Montaldi and Guttmann and Prellberg, we can write for 
the d-dimensional hyper-cubic LGF 

[ww-^ 2 ;) £ (iT^) 2 <») 

V 7 k 1 +k 2 + ...k d =n V a 7 

The 5 dimensional LGF satisfies a 5th order ODE of degree 3, being equation number 
188 on the list [4j, and which can be "pulled back" to a degree 12 ODE of 4th order, 
that is also C-Y. 

For d < 5 the multinomial expression can be simplified to products over binomial 
coefficients. Specifically, writing 

P($*) = £«*(^) n > (36) 
we have for d = 2, 3, 4 respectively 
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No such simplification is known for d > 5. 
However in [2 9) I give the recurrence 



2 

n ' 



w = e ; ) ^ 

where 

S { n } = E 



k 1 lh 2 l...k d l, 

So for d — 5 we can write the multinomial as a double sum: 



si 5) = E E 

fc 1= fc 2 =0 



n\ \ (2ki \ (2k 2 



k\\k 2 \(n — k\ — k 2 )\ / \ k\ J \ h 



>2 



Hence 



t") E (/,,1/J.'.. /,-,;) ^ 

fci+...fc 5 =n 

n J ^ 1" \kx\k 2 \{n -kt- k 2 )\) \ki)\ha 

fcl— U K2— U 

Joyce [35 j has also studied the LGF of the hypercubic lattice, and obtained a rapidly 
convergent asymptotic expansion for the coefficients of the d-dimensional lattice LGF, 
as well as elucidating the behaviour of the LGF in the vicinity of the dominant branch- 
point singularities. 

Finally we note the identity for the following integral of the d th power of the mod- 
ified Besscl function Iq. 



r e -^ { 2 Z t)dt=±f^) £ 



k 1 J rk 2 J \-...k d — n 



k x \k 2 \...k d \ 



2 

2n 



(38) 



This is just the generating function for returns on the d-dimensional hyper-sc lattice. 
The corresponding integral with Iq replaced by Jo produces the same result with an 
extra factor (—1)™ in the sum on the r.h.s. 

2.3 Hyper diamond lattice 

There is some confusion in the literature as to the 4d generalisation of the diamond 
structural. 

There are two types of site, and the relevant unit vectors are: 
(±1,0,0,0), (0,1,0,0), (0,0,1,0), (0,0,0,1) 



6 I too was confused, and my confusion was resolved by David Broadhurst. A correct discussion can 
found at www.hermetic.ch/compsci/lattgeom.htm 
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and 

(±1,0,0,0), (0,-1,0,0), (0,0,-1,0), (0,0,0,-1) 

Calculating the structure function by taking the Fourier transform of the individual 
step probabilities gives 



A 2 = 2 cos(fc 2 — fc 3 )+2cos(A:2— k 4 )+2 cos(fc 3 — fc4)+4cosfci(cos &;i+cos&2+cos/c3+cos A: 4 )+3. 

To evaluate the LGF, I did one of the integrals exactly, and generated the series 
expansion by doing the remaining 3-dimensional integral within Maple. This is far 
from the most efficient method, but was adequate for our purpose, as the series for 
the number of 2n step returns was immediately recognizable as the coefficients of the 
square of the 5-dimensional multinomial coefficients. 



i+j+k+l+m=n v J 



2 

2n 



As the corresponding result for the 3-dimensional diamond lattice, and its two- 
dimensional counterpart, the hexagonal lattice, are given by the squares of the four- 
dimensional and three-dimensional multinomial coefficients respectively, this is hardly 
surprising. 

I gave the ODE for this in a 1993 paper with Thomas Prellberg [29], and it is simply 
related to the gen. fn. for 5d staircase polygons! It is T>P(0, 5z) = where 

V = 8 4 -z(356 i +706 3 +636 2 +286+5)+z 2 (0+l) 2 (2596 2 +5180+285)-225z 3 {6+l) 2 (9+2) 

and is number 34 on the list |3] of 4th order, third degree Calabi-Yau equations. 
So using the results of the previous subsection, we immediately have 



v ( n - ) 2 = v "v ( nl \ 2 ( 2kl \ ( 2k2 



(39) 

Higher dimensional generalisations are obvious. The <i-dimcnsional hyper diamond 
LGF has coefficients given by the sum of the squares of the (d+ l)-dimensional multi- 
nomial coefficients, 



,(<*) = 

ki+k.2+...k d+1 - 



■ - (wsrrwJ- (40) 



The 5th order differential equation for the d = 5 case is given in [29]. It is also given 
as ODE number 130 in the list [4], where it is shown that it can be "pulled back" to a 
4th order CY ODE of degree 12. 

Broadhurst [13] has also calculated the Yukawa coupling for the LGF ODEs of the 
diamond lattices in dimensions 3 < d < 10 and extracted the instanton numbers nk{d) 
of the Lambert series. By numerical experimentation he found that n,k(d)/k 2 is an 
integer for d > 3, k > 0, and gave explicit polynomial representations in d for rik(d) for 
k < 7, and verified the conjecture that these are polynomials of degree k up to d = 9. 
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The (i-dimensional diamond lattice Green function can also be expressed compactly 
in terms of Bessel function integrals, as shown by Glasser and Montaldi [25] • In fact 



/•OO °° 

/ tI i d+1) (zt)K (t)dt = P d (d-(d+l)z) = Y / E 

•* „_ n u. ■ j — ■ h, 



>i!fc 2 ! . • ■ fcd+i!. 



n=0 fci+fe 2 + ...fc d +l = 



2 



(41) 



2.4 Connection between the hyper diamond and hyper cubic 
lattices 



We have seen that the (i-dimensional hyper-diamond lattice LGF has coefficients (|40| 
given by the sum of the squares of the (d + l)-dimensional multinomial coefficients. 
The (i-dimensional hyper cubic lattice has LGF coefficients ([33)) given by i 2 ™) times the 
sum of the squares of the d-dimcnsional multinomial coefficients. Thus there appears 
to be a simple relationship between the LGFs of the two lattices. 

As it happens, this has been formalised in another context by Guttmann and Prell- 
berg [59] and Glasser and Montaldi 25 . In [21] a study was made of the generating 
functions whose coefficients were the squares of the d-dimensional multinomial coeffi- 
cients, as part of a study of d-dimensional staircase polygons. For d — 4 this generating 
function was shown to be connected to the LGF of both the (three dimensional) f.c.c. 
and diamond lattices. 

In [25] this was clarified, as Glasser and Montaldi showed that the generating func- 
tion for the (i-dimensional multinomial coefficients, 

U^)= E (\ + ---\ kd \ 2 x ^+-^ (42) 
k u ...,k d =0 V fel '"-> fc < i / 

is equivalent to the generating function for the (i-dimensional hyper cubic LGF 

p r \ _ J_ f°° [°° de 1 ---de d 

d[Z> Wo "Jo l-|[cos(^) + ---+cos(^)] [6> 
through an Abel transform. In particular, they showed that 

* Z d (t 2 z 2 /d 2 ) 



Wo VI - 1 2 



(44) 



and 

rl tP d (dxt) 



dt . (45) 

Here we capitalise on this result by pointing out that Z d+ i(x 2 ) is the LGF for the 
(i-dimensional hyper diamond lattice, with coefficients given by eqn (|40p . 

Invoking this connection at the coefficient level, implementing eqn (|44p only requires 
the calculation of the integral 

t 2n , 1 (2n 
dt = — 
VT^i2 4™ 
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which reveals the origin of the prefactor ( 2 ^) in the expression for the coefficients of 
the hyper cubic LGF. 

Finally, comparison of eqns (|4Tl [38l l44l) yields the Bessel function relation 

r\ p 1 J |>00 

e- t I${zt/d)&t= - / tI$(ztu/d)K (t)dt. (46) 

t Jo Vl-u 2 Jo 

Other Bessel function identities discovered through their connection with LGFs can be 
found in Broadhurst [15] . 



2.5 Hyper face-centred cubic lattice 

The four-dimensional hyper-fcc lattice Green function is given by 

„,* , 1 r r dkidk 2 dk 3 dk A 
P(0;z) 



(t) 4 Jo Jo 1 - f A 

where A = {c\Ci + c\c-$ + C1C4 + C2C3 + C2C4 + C3C4) and q = cos/c^. 

In [27 I found the CY ODE satisfied by P(0; 24z) by doing two of the four inte- 
grations exactly, and then expanding the integrand in a power series and integrating 
term-by-term within Maple. In this way I identified the ODE from a series of some 41 
terms, which contains barely any confirmatory coefficients. Subsequently Broadhurst 
[14] obtained the series to 100 digits, thus providing ample confirmation of the result. 

The two integrals were performed as follows: 

Set a = 1 — |(c 2 c 3 + c 2 C4 + C3C4); b = |(c 2 + c 3 + c 4 ). 

Then the integrand is [a — frcosfci] -1 . 

Use 

1 r de 1 1 



a - b cos 6 yja 2 - b 2 y/(a + b)(a - b) 

to eliminate k\. Next write (a + b)(a — b) = e(c — cos k2){d — cos ^2), where c, d, e are 
independent of ki , and use 

d6 2K(fc) 



y/ {c - cos e)(d- cose) i/(c- l)(d+l) 

to eliminate fc2, where fc 2 = rg~w^x] ■ We are left with a two-dimensional integral, 
which was expanded as a power series in z and integrated term-by-term in Maple. We 
got to 41 terms in a few hours, then searched for an ODE. With 6 = z-4-, the LGF 
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satisfies £>P(0; 24z) = 0, where^ 

V = 9 i + z(396» 4 - 306> 3 - 198 2 - 4(9) 

+ 2z 2 (16<9 4 - 10706» 3 - 10576 2 - 6760 - 192) 

- 36z 3 (1716» 3 + 5666> 2 + 6006> + 316) (30 + 2) 

- 2 5 3V(+3846> 4 + 15426» 3 + 26356> 2 + 21730 + 702) 

- 2 6 3 3 z 5 (13936» 3 + 55716> 2 + 83786* + 4584) (1 + 6) 

- 2 10 3 5 z 6 (316» 2 + 1056> + 98) (1 + 6)(6 + 2) 

- 2 12 3V(6» + l)(6> + 2) 2 (6 , + 3). 

This is a 4th order, degree 7 Calabi-Yau ODE with regular singular points at 0, 1/24, 
-1/3, -1/4, -1/8, -1/12, and oo. It was new at the time it was found, and is one of only 
a few known C-Y ODEs of degree 7. It is now ODE number 366 in the list of Almkvist 
et al. g]. 

We do not yet have a nice expression for the series coefficients in terms of binomial 
coefficients. A result is given in [3J, which we have simplified a little, giving 

(2i\ /2A [2k\ (I + m\ (2(1 + m 

ii \ I n — 21 — 2m \ / 2i — I — m 
2{l + m))\n-2i-l-m)\ i-k-l 

In two dimensions there is a simple connection, at the coefficient level, between the 
coefficients of the triangular lattice LGF and the honeycomb lattice LGF, given by eqn 
(fTTj) . In three dimensions a similar connection exists between the coefficients of the 
fee and diamond lattices, see eqn (|17|) . Unfortunately there seems no such connection 
between 4d fee and 4d diamond lattices. 

Subsequently David Broadhurst pointed out to me that a computationaly far more 
efficient way to obtain the series is to expand the integrand directly. The coefficients 
of z n are given by readily evaluatable multinomial coefficients, and the integrals are 
then just integrals of powers of cosines over a half-period, and we know 

i^°°cos 2 "xdx= ( 2 n n )/4", 

so the general term can be written down as a five-fold sum over products of a multi- 
nomial coefficient and 6 binomial coefficients. He also showed [14] that the degree 
of the ODE could be reduced to 6 by considering the auxiliary function Fa,{z) = 
F±[z/(1 — 18z)]/(l — 18z) where F^[z) is the 4-dimensional hyper fee LGF. The re- 
sulting ODE has the C-Y properties. In 14 Broadhurst also calculated the first 100 
instanton numbers. The first few of these are 3, —4, 64, —253, 4292, —25608. 

Indeed, Broadhurst in a heroic calculation [2] also found the ODE for the 5- 
dimensional hyper fee LGF. Disappointingly, it destroys the regularity seen to date, as 



7 This corrects a sign error in the coefficient of z in [27] . Also, in the list of singularities given there, 
-1/18 should be replaced by —1/3. 
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it is not a 5th order but a 6th order ODE, of degree 13, and lacks the MUM property. 
By a similar transformation to that used for the 4-d fee LGF, the degree can be reduced 
to 12. Attempts at factorizing the differential operator have been unsuccessful. 

2.6 Lattice Green function type integrals 

There are several other structure functions that are "obvious" generalisations of the 
3d lattice structure functions, and these have been investigated similarly. 
The first is: 

1 r r dh dk 2 dk 3 dk A ^ 

mz)= wl -I i-zx ^; nZ ' 

where A = (C1C2C3C4 + S1S2S3S4) and Cj = cosfci, Sj = sinfej, Then I find that P(0; z) 
is exactly the same as the LGF of the 4-d hypercubic lattice eqn flMI) . with coefficients 
a n given by: 




A more interesting structure function is A = (C1C2C3 + C1C2C4 + C1C3C4 + C2C3C4). 
Again we can efficiently calculate the Green function (I drop the adjective "Lattice" as 
there is no obvious lattice driving this structure function) by expanding the integrand 
and integrating term by term. This yields a four-fold sum of products of a multinomial 
coefficient and four binomial coefficients. Broadhurst, and independently I, showed that 
the resulting series satisfies an ODE of order 8 and degree 16. The regular singular point 
at the origin has roots of the indicial equation equal to (4 times), |, |, |, (2 times), 
so is not MUM. As this structure function does not correspond to any known lattice, 
that absence of the MUM property is of no concern. 

3 Constant term (CT) identities 

An alternative way to view the series expansions for the number of distinct returns 
to the origin of a random walker after n steps, which is given by the (appropriately 
scaled) LGF, is to ask for a constant term formulation. The CT is manifestly a more 
obviously combinatorial approach than the integral representation of the LGF. The CT 
enumerates the number of distinct n-step returns to the origin by directly generating 
all possible sums of lattice basis vectors that sum to zero. It is also immediate that the 
value is an integer - something that is not so obvious from the integral formulation. 

The idea here is to write down a function which when raised to the n th power has 
a constant term equal to the number of n-step returns. For the d-dimensional hyper 
cubic and hyper b.c.c. lattices, only even values of n are used, as returns to the origin 
must have an even number of steps. For the triangular and hyper f.c.c lattices, all 
values of n are used. For lattices with two types of site, such as the honeycomb and 
diamond lattices, all values of n are used, but the result is the number of 2n step 
returns to the origin. 
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For example, for the square lattice the required expression is 

f(x,y) = (x+-)(y + -). 

x y 



Then CT f(x,y) 2n = ( 2n ) 2 . Alternatively, if 



/ \ , 1 Is 

g{x,y) = (x-\ hy+-j, 

x y 

thenCT 5 (x,y) 2 " = ( 2 ;) 2 . 

For the triangular lattice the relevant function is 

, . . 1 1 y x. 

h(x, y) = (x+- + y+- + - + -). 

x y x y 

For the honeycomb lattice one has 

u(x, y) = {l + x + y)(l + 1/x + 1/y) 
For the diamond, s.c, b.c.c, f.c.c. lattices respectively the relevant functions are: 
(l/x + x + z{y + l/y)){x + 1/x + (y + 1/y)/ z) diamond 

.1 1 1. 

[x-\ \-y-\ hz + - s.c. 

x y z 

(x + -)(y + -){z + -) b.c.c. 

x y z 

(x + -)(y + -) + + -){z + -) + (z + -){y + -) f.c.c. 

x y x z z y 

Some of these results were first given by Domb in [20] . 

For the four dimensional diamond, s.c, b.c.c, f.c.c. lattices respectively the relevant 
functions are: 

(1/x + x + zy + z/y + w/x)(x + 1/x + y/z + 1/yz + x/w) 4d diamond 

1 1 1 In aj 

(x H hyH hzH hioH ) 4d s.c. 

a: y z w 

(a;+ -)(y + -)(z+ -)(w + — ) 4d b.c.c. 

a; y z w 

(x+-)(y+-) + (x+-)(z + -) + (z + -)(y+-) + (w + —){x+-+y+- + z + -) 4d f.c.c. 

x y x z z y w x y z 

These then give constant term formulae for the C-Y ODEs numbers 34, 16, 3 and 
366 respectively in the list of Almkvist et al. [3]. 
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3.1 CT formulations and Mahler measure 



There is a significant, but not yet fully understood connection between the CT formu- 
lation of LGFs and the Mahler measure of a polynomial. Papers developing aspects 
of this connection have been brought to my attention by Almkvist, who mentioned 
Stienstra's work [321133], Glasser who mentioned Boyd's paper [T2] and Zucker who 
drew my attention to the recent paper by Rogers |41j . 
If one writes a polynomial as 

p(z) = a(z - a\){z - a 2 ) ■ ■ ■ (z — a n ) 

then the Mahler measure M(p) of the polynomial is given by the product of the absolute 
value of a and those roots with magnitude at least 1. That is to say, 

M(p) = \a\ J] |a,|. 

|a;|>l 

More generally, for a Laurent polynomial F(x,y,z) one normally defines the loga- 
rithmic Mahler measure m(_F), which we denote Mm, and the Mahler measure M(_F), 
respectively as 

m ( F ) = 7^l/// \og\F(x,y,z)\ — ^- — (48) 

(2ttz) 3 / / J lxHyHzl=1 x y z 

and 

M(F) = exp(m(F)), (49) 

where the generalisation to more variables is obvious. 

It is not our purpose here to summarise the contents of the cited papers, but we 
mention those of direct relevance to this work. Boyd [H] considers the IMm of the 
Laurent polynomial 

(aH h y H hi 

x y 

which is, apart from the additive factor 1, the CT kernel for the square lattice LGF 
given above. Boyd investigates generalisations of the conjectured result of Deninger 
[19] that 

1 - — — 

y + ' ~ (2tt)2 

where L(E, s) is the L-function of the elliptic curve E of conductor 15 that is the 
projective closure of (x + - + y + - + 1). (In fact Dcninger's conjecture included a 

x y 

constant multiple, which Boyd has shown to be precisely 1 to more than 25 digits). 
Boyd goes on to conjecture a number of other evaluations of the IMm of two-variable 
Laurent polynomials in terms of L'(E,0), where E is a given elliptic curve. 
Rogers [41] considers two three- variable Laurent polynomials, and studies 

1 1 1 

h(u) := mix H hyH h z -\ \-u), 

x y z 



m(x+-+y + z + l) = j^L(E, 2) = L'(E, 0), 
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and 



g(u) := m (4 - u+ (y + -){z + -) + (a: + -)(y + -) + (x + -){z + -)) , 
\ y z x y x z J 

which are clearly related to the CT kernel of the simple cubic and face-centred cubic 
lattices respectively. By considering these and other Laurent polynomials, and their 
Mm, Rogers establishes a number of new identities for the hypergeometric function 
5-F4, and new formulae for l/n, some of which we have mentioned above in Section 1.3. 

Other results of Rogers can be transferred to our problem. For example, with a 
trivial change of variable, Theorem 3.1 of Rogers can be written 

1 /l 1 2 27z 4 \ 

P(0,z)diamond = {l _ z 2/ A f F ^ [ 3 ' 2 ' 3 5 ^ 64(1 - Z^/Af J ' ^ 

Given the known relationship between the diamond and f.c.c. LGFs, we can also write 
P(0, z) fcc = 3 F 2 Q, i, 1 1, 1; f!fcfl j . (51) 

Alternative expressions in terms of the squares of a 2F1 hypergeometric function 
were first given by Joyce in 1994 [32], who also gave expressions in terms of different 
3F2 hypergeometric functions. This would give a slightly different version of Theorem 
3.1 of Rogers. The expressions above for the LGFs on the f.c.c. and diamond lattices 
are slightly simpler than those in [32]. 

Finally, Stienstra's two rich papers [52] 03J make connections between Mahler mea- 
sures, the partition function of some dimer models, and instanton numbers in string 
theory. It is known that the square lattice dimer partition function can be related to 
a certain class of spanning tree generating function, which can in turn be obtained by 
integration of the LGF [28] . We will expand on these connections in a future publica- 
tion. 



4 Conclusion 

We have summarised known results for LGFs in lower dimension, and present less- 
well-known results for higher dimensional results. These are then combined to give a 
cohesive view for some lattices as a function of dimensionality (for the <i-dimensional 
diamond, simple-cubic and body-centred cubic lattices). For the face-centred cubic 
lattice, the situation is less satisfactory, though results for 2 < d < 5 are given, but we 
cannot give an expression that holds for arbitrary dimensionality, except through the 
constant-term formulation. 

For two dimensional lattices, the LGF can be expressed in terms of the complete 
elliptic integral K(fc). For three dimensional lattices, the LGF can be expressed in 
terms of products of two complete elliptic integral K(fc_)K(fc + ), or as the square of 
a 2-F1 hypergeometric function or as the first power of a 3-F2 hypergeometric function. 
We find slightly simpler 3F2 representations for the diamond and f.c.c. lattices. 
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Joyce [35] also found a remarkable result in that it is possible to express the LGF 
for the four common three-dimensional lattices discussed above in terms of the LGF 
for the two-dimensional honeycomb lattice, ([7]). First, define 

R(0= Phoney @,3z). 

The series expansion of has integral coefficients giving the number of 2n step 

returns to the origin of a random walker on the honeycomb lattice. Then 

Pfcc&z) = p/ cc (o,o = (i - sennm 2 

where 

z = z(£) = -12£ 2 /(l-3e 2 ) 2 - 
For the simple-cubic lattice Joyce finds 

P sc (0,z)EEP sc (0,O = (l-9e 4 )[i?(O] 2 

where 

z 2 ee z 2 (£) = 36£ 2 (1 - 9£ 2 )(1 - £ 2 )/(l - 9C 4 ) 2 - 
For the b.c.c. lattice he finds 

p bcc (o,z) ee A cc (o,o = (i - 9e 2 ) 1/2 (i ~e) 3/2 [R(o? 

where 

z 2 ^z 2 (0 = -64£ 6 /[(i-9£ 2 )(l-a 3 ]- 
Finally for the diamond lattice Joyce gives 

p d (o, z) ee p d (o,o = (i - 9^ 2 )(i - emm 2 

where 

« 2 = z a (0 = -i6e a /[(i-9e a )(i-e a )]- 

For four dimensional lattices, the LGF can be expressed as an integral with an 
integrand that is a product of an algebraic function and a product of two complete 
elliptic integrals K(fe) for the s.c. lattice by combining eqns (fT2j) and (|44|) . giving 



4d 



1 K(fc+)K(/fc_)dx 



2 



a similar representation can be found for the 4d b.c.c. lattice. It remains to be seen 
whether this can be achieved for the 4d f.c.c. and diamond lattices. More generally, 
it still remains to be seen if the results for the 4d LGFs can be presented in a unified 
manner, as Joyce has done for the 3d lattices, as described in the preceding paragraph. 
What is clear however is that the remarkable results that pertain for the 3-d lattices 
are due to the fact that the underlying 3rd order ODE has the almost-magical Appell 
reduction property, allowing the solutions to be expressed in terms of the solutions 
of an associated second-order ODE. For the 4-d LGFs, we have the almost equally 
remarkable property that the underlying ODEs are all of Calabi-Yau type. For the 5-d 
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LGFs, this beautiful property appears to be broken by the 5-d f.c.c lattice, though it 
holds for the other three lattices, all of which have solutions that can be expressed in 
terms of an associated fourth order ODE. 

For the d-dimensional s.c. lattice, we can write the result for arbitrary dimension- 
ality in terms of an integral over the d th power of a modified Bessel function, given by 
(|38|) , with a corresponding result for the e?-dimensional diamond lattice, given by (|41[) . 
The corresponding result for the ci-dimensional b.c.c. lattice is even simpler, being 
expressible in terms of a dFd—i hypergeometric function, or more simply still we can 

write the number of 2n step returns as ( 2 ^) . We have no corresponding result for the 
hyper-fcc lattice. 

For the two-dimensional triangular lattice, the coefficients can be expressed in terms 
of those of the honeycomb lattice, as shown in eqn (fTTj). For the three-dimensional f.c.c. 
lattice, the coefficients can be expressed in terms of those of the diamond lattice, as 
shown in eqn (fT7|) . Unfortunately, it seems that for the 4d f.c.c. lattice, there is no 
simple expression in terms of the coefficients of the 4d diamond lattice. More precisely, 
we have been unable to find one. 

The number of returns has been given by a constant-term formula for all lattices, 
incidentally augmenting the table in [|] , and it is only this formulation that allows us 
to generalise the f.c.c. lattice to arbitrary dimension. That is to say, we can write 
down an expression for the number of returns of arbitrary length on any regular lattice 
of ci-dimensions. Unfortunately, this still leaves one with a non-trivial computational 
problem in seeking the numerical values. 

We also find that some remarkable Ramanujan-type formulae for 1/n involve the 
coefficients of the 3-dimensional diamond, s.c and b.c.c LGFs. 
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Appendix A— Yifan Yang pullback 

In this appendix we discuss the notion of YY-pullback (a contraction of the eponymous 
Yifan Yang pullback) . In preparing this appendix we have drawn heavily on the papers 
by Almkvist [lj and by Almkvist, Van Straten and Zudilin [5]. First, recall the result 
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of Appell, [7], established 130 years ago, that if a 3rd order ODE can be written as 

/"'Or) + 3P(x)f"(x) + [2P 2 (x) + P'(x) + 4Q(x)]f'(x) + [4P(x)Q(x) + 2Q'(x)]f(x) = 0, 

(52) 

then its solution can be expressed in terms of the two linearly independent solutions, 
gi{x) and g 2 {x) of the 2nd order ODE 

g"(x) + P(x)g'(x) + Q{x)g(x) = 0. (53) 

In fact 

f(x) = Agl(x) + B 9l (x)g 2 (x) + Cg\{x). 

Equation (|52l) is called the symmetric square of eqn (|53l) . An example of the utility 
of this connection can be found in .29 , while another example establishes Clausen's 
formula relating the square of a certain 2-F1 hypergeometric function to a 3-F2 hyper- 
geometric function. 

Consider now a 4th order ODE satisfying the C-Y conditions described by eqn (|3"2"j) 
and the four other conditions given in Section 2. 

h {m \x) + P(x)h"'(x) + Q(x)h"(x) + R(x)h'(x) + S(x)h{x) = 0. (54) 

The CY condition implies that the four linearly independent solutions /io, hi, h 2 , /13 
of (|54"|) are as in (|31l) . Then the six Wronskian functions 

Wjk = W(hj ,h k )= det *£),0<j<k< 3, (55) 

are not linearly independent over C. Indeed, 

W03 = W12, 

a condition equivalent to eqn (|3"2"1) . The Wronskians xwij satisfy a 5th order ODE: 

u {v) (x) + P{x)u (lv) {x) + Q(x)u"'(x) + R(x)u"(x) + S(x)u'(x) + f{x)u{x) = 0. (56) 

Let wq = xwqi and wi = xwq2- Then if we construct the Wronskian of these Wron- 
skians, we have 

W(w Q ,wi) = x 2 ?/oexp(-i J P(x)dx). 

This allows one to recover the solution ho from the solutions wq, w\ of the fifth order 
equation by using 

2 2 - 
P(x) = - + -P(x). 

x 5 

Then multiplying by a suitable factor we get the Yifan-Yang pullback: 

y(x) = x 5/2 y / W(w ,w 1 )exp(-^ J P(x)dx). 

The motivation for this is given in [2], but the point is that the YY-pullback is 
usually of lower degree (typically half the degree) of the ordinary pullback found in the 
table @]. 

A number of examples are given in [3] and a few are cited in the text above. For 
example, the 5th order ODE satisfied by the 5-dimensional hypercubic LGF is of degree 
3, as shown in entry number 188 of [4] but its fourth order YY-pullback is of degree 6, 
though is not yet in the Table 1 pQ. 
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Appendix B— Calabi-Yau differential equations 



In this appendix we give an attenuated account of the much richer and more detailed 
material by Almkvist, van Enckevort, van Straten and Zudilin in jH 0] and [5]. The 
name Calabi- Yau differential equations comes from a search by the just-cited authors 
for "arithmetically nice" differential operators of order 4 and 5. Before discussing what 
is meant by the adjective, first, following [5], consider the second order, second degree 
linear ODE, 

{9 2 - z{a9 2 + a6 + b) + cz 2 (6 + l) 2 )f(z) = (57) 

where, as usual, 9 — z4-. We require it to have a unique analytic solution fo{z) = 
1 + a n zTl where the coefficients a n are integers. There will also be a second solution 
fi(z) = fo(z) logz+<?(z), where g(z) is an analytic function and g(0) = 0. Furthermore, 
we require the function 

q(z) = zcxp(g(z)/ f (z)) = b n z n 

n>l 

to also have coefficients b n that are integers. In Zagier [Uj, the result of an exhaustive 
search for values of {a, 6, c} satisfying these conditions is reported. A necessary condi- 
tion is that the values {a, 6, c} are integers. However there are only 14 non-degenerate 
sets of triples found. 

Moving to third order ODEs, if we consider the third order, second degree linear 
ODE, 

{9 3 ~ z(29 + \){a9 2 + a9 + b) + cz 2 {9 + l) 3 )/(z) = 0, (58) 

and if we impose exactly the same conditions as for the second-order equation, it was 
found by Almkvist et al. [4] that there are again only 14 triples {a, 6, c} satisfying 
these conditions. As an aside, it is pointed out in [6] that the second-order equation 
with the triple {11, 3, —1} arises in Apery's proof of the irrationality of C(2), while the 
third-order equation with the triple {17, 5, 1} arises in Apery's proof of the irrationality 
of C(3). Furthermore Almkvist et al in [5] prove a natural bijection between the set 
of triples {a, 6, c} and the set of triples {a, 6, c}, and use that observation to link the 
solutions by a rational map. 

In order to generalise these results to higher order ODEs (specifically orders 4 and 
5), the approach taken in [6l|4] was to require the differential operator to satisfy simi- 
lar conditions to those satisfied by the second- and third-order ODEs just considered. 
These are: 

(i) the ODE must be Fuchsian (all singular points are regular), 

(ii) it has the MUM property (all exponents are zero for the singular point at the 
origin), 

(iii) the unique analytic solution at the origin, fo(z) = 1 + ^a n z n has integral coeffi- 
cients a„, 

(iv) the solution f\{z) — fo(z) logz + g(z), g(0) = 0, gives an expansion with integral 
coefficients b n for the function 

q(z) = zexp(g(z)/f (z)) = ^ b n z n . 

n>l 
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The Calabi-Yau condition (|32p was additionally imposed in [51 0] , though in hind- 
sight it appears, at least experimentally, to be a consequence of the four conditions 
above. Likewise the integrality of the (possibly scaled) instanton number^. By a 
variety of ingenious constructions and exhaustive searches, there are now more than 
350 known equations of degree 4 catalogued in g]. They are called Calabi-Yau since 
they can often be identified with Picard-Fuchs differential equations for the periods of 
1-parameter families of Calabi-Yau manifolds. It is believed that all these ODEs are of 
geometric origin. Technically, and here we quote [5], this means that they correspond 
(as subquotients of the local systems) to factors of Picard-Fuchs equations satisfied by 
period integrals for some family of varieties over the projective line. 

Taking this connection with geometry a little further, Almkvist et al. show that, 
of the 14 solutions in the case of second order ODEs, four are hypergeometric, four 
are Legendrian and six are sporadic. The hypergeometric solutions correspond to the 
value c = in ([ST]) . In the Legendrian case one has c = a 2 /4 in ([ST]) . In both these 
cases the differential operators are Picard-Fuchs operators of the extremal rational 
elliptic surfaces with three singular fibres [35]. Finally, the six sporadic cases are also 
of geometric origin, arising as Picard-Fuchs equations of the six families of elliptic 
curves with four reduced singular fibres [9] [38] . 
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